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Abstract. The problem of backward dynamics over the ring of p-adic 
integers is studied. It is shown that Inverse Limit Theory provides the 
right framework. Backward iterations of a polynomial with p-adic integer 
coefficients are constructed by solving congruences modulo powers of p, 
which inturn are solved by Hensel's lifting lemma. 
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1 Introduction 

Dynamical systems originally arose in the study of systems of dif- 
ferential equations used to model physical phenomena. One simpli- 
fication in this study is to discretize time, so that the state of the 
system is observed only at discrete steps of time. This leads to the 
study of the iterates of a single transformation. One is interested 
in both quantitative behavior, such as the average time spent in a 
certain region, and also qualitative behavior, such as whether a state 
eventually becomes periodic or tends to infinity. 
A discrete-time dynamical system consists of a non-empty set X and 
a map f : X ^ X. For n G N, the nth iterate of (p is the n-fold com- 
position ^t"! = 0000.. .00 (n times); 0'°] is defined to be the identity 
map. If is invertible, then 0"" = 0~^o0~^o...o0~^ (n times). Since 
0[n+m] _ 0[n] ^ thesc iterates form a group if is invertible, and 
a semigroup otherwise. For a given a in X, the forward orbit of a is 
the set 

O^ia) = {0"(«) : n > 0} 

If the orbit C</,(a) is finite then a is said to be a pre-periodic point, 
otherwise a is said to be a wandering point. 

The central problem in dynamics is to classify the points a in the 
set X according to the behavior of their orbits 0^(^a)- In practice X 



2 Backward Dynamics and p-adic Approximations 



usually has additional structure that is preserved by the map /. For 
example, (X, /) could be a measure space and a measure preserving 
map; a topological space and a continuous map; a metric space and 
an isometry; or a smooth manifold and a differentiable map; a finite 
set (e.g., Finite Field) and a polynomial. 

1.1 Symbolic dynamics 

Symbolic dynamics arose as an attempt to study such systems by 
means of discretizing space as well as time. The basic idea is to di- 
vide up the set of possible states into a finite number of pieces. Each 
piece is associated with a "symbol" , and in this way the evolution of 
the system is described by an infinite sequence of symbols. This leads 
to a "symbolic" dynamical system that mirrors and helps us to un- 
derstand the dynamical behavior of the original system. Computer 
simulations of continuous systems necessarily involve a discretiza- 
tion of space, and results of symbolic dynamics help us understand 
how well, or how badly, the simulation may mimic the original. Sym- 
bolic dynamics by itself has proved a bottomless source of beautiful 
mathematics and intriguing questions. As polygons and curves are 
to geometry shift spaces arc to symbolic dynamics. The set 

E ^ {0,1,2,..., 

is called the sequence space on m symbols 0, 1, m — 1. 
The most important ingredient in the sequence space is the shift 
map (7. The shift map a : E ^ E is given by cr((so, Si, S2, ■■)) — 
(•Si, ^2, -Sa, •••)• 

The shift map discards the first entry in the sequence and shifts all 
other entries one place to the left. 

The distance between two sequences s — {sq, Si, ...) and t — {to,ti, ...) 
is given by rf(s,t) = E::o^- 

For an integer m > 1, set = {1,2, ...,m}. Let = be the 
set of infinite two sided sequences of symbols in and Z'^ = A^ 
be the set of infinite one-sided sequences. The pair {E^, (j) is called 
the full two sided shift; {E:^, a) is called the full one sided shift. The 
two-sided shift is invertible. For a one-sided sequence, the leftmost 
symbol disappears, so the one-sided shift is non-invertible, and ev- 
ery point has m pre-images. Both shifts have m" periodic points 
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of period n. The shift spaces are compact topological spaces in the 
product topology. This topology has a basis consisting of cylinders 

^jT,''jr={^ = (^') • = -^i'^ = 1,2, rii < n2 < ... < Uk 

are indices in Z or N, and ji G A^. Since the preimage of a cylinder 
is a cylinder, a is continuous on and homeomorphism on Em- 
The metric d{x,x') = 2~\ for / = min{\i\ : Xi 7^ x'^} generates the 
product topology on 1^+ and Em- In the product topology periodic 
points are dense and hence there are dense orbits. A detailed study 
can be found in [1], [5]. 

2 Arithmetic Dynamical Systems 

Classically, discrete dynamics refers to the study of the iteration of 
self-maps of the complex plane or real line. Arithmetic Dynamics 
is the study of number theoretic properties of dynamical systems. 
Arithmetic dynamics is discrete-time dynamics (function iteration) 
over arithmetical sets, such as algebraic number rings and fields, 
finite fields, p-adic fields, polynomial rings, algebraic curves, etc. A 
thorough introduction is given in [1], [3]. 

In this paper we study the problem of backward dynamics of any 
polynomial of finite degree over finite rings Z/p"Z using Hensel's 
lifting lemma, as usual p denotes a prime. 
Backward-iteration sequences given by 

are of a different nature because a point could have infinitely many 
pre-images as well as none. If the given forward moving map is a 
quadratic map, the corresponding backward map is a square root 
map; if the given map is a cubic map the corresponding backward 
map is the cube root map and so on. Thus essentially we solve 
f{x) = as a polynomial over the defining set. For e.g., the Julia set 
can be found as the set of limit points of the set of pre-images of (es- 
sentially) any given point. Unfortunately, as the number of iterated 
pre-images grows exponentially, this is not feasible computationally 
when the underlying set is the set of Real or Complex numbers. 
In general, maps of higher degree (> 5) are not suitable for backward 
dynamics over M or C. As there is no explicit formula for solving a 
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polynomial of degree > 5, roots can be found by using standard tech- 
niques from Nummerical Methods over M or C and we either arrive at 
a null sequence or constant sequence after some backward iteartion. 
In such cases nothing can be said about the behavior of trajectories. 
But when the set is finite (endowed with algebraic structure) and the 
map is a polynomial it is possible to retrieve the pre-images (roots), 
if they exist, with respect to different prime power moduli and thus 
study the structure of pre-images locally at that prime. 
We show that this problem can be well understood over the ring of 
p-adic integers Zp. This process deals with an imporatnt branch of 
mathematics called Inverse Limit Theory. Below we discuss some 
basics of Inverse Limit Theory and the p-adic integers. 

2.1 Inverse Limits 

Let Xo,Xi,X2,... be a countable collection of spaces, and suppose 
that, for each n > 0, there is a continuous mapping : X„ — )■ 
The seqence of spaces and mappings /„} is called an inverse 
limit sequence and may be represented as 

fn + l -tA fn -1^ /2 -lA /l -lA 

... > J^n ^ — > — > 

Clearly, if n > m, there is a continuous mapping : Xn — )■ Xm 
given by the composition fn,m = fm+i ■ fm+2 ■ ■ ■ fn-i ■ fn- 
Consider a sequence {xq,Xi, ..) such that each point of 

the space X„ and such that Xn = fn+i{xn+i), n > 0. Such a sequence 
can be identified in the product space H^^o"^" considering a 
function (p from the nonnegative integers into H^o"^"' given by 
(j){n) = Xn- Thus the set of all sequences is a subset of H^o"^'! 
has a topology as a subspace. This topological space is the inverse 
limit space of of the sequence /„} denoted by X = ^m(X„, /„). 

Theorem 1. Suppose that each space Xn in the inverse limit se- 
quence {Xn, fn} is a compact Hausdorjf space. Then X is not empty 
If. 

Theorem 2. A space X is a compact Hausdorff space with dim{X) < 
if and only if X is an inverse limit of finite discrete spaces 
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A finite discrete space is totally disconnected, compact and Hausdorff 
and all those properties carry over to inverse limits too. 
A detailed study of Inverse limit spaces can be found in [B] . 

2.2 p-adic Integers 

An important class of such inverse limits is given by rings of p-adic 
integers Zp. For every n > 1, let An = Z/p"Z. An element of An 
defines in an obvious way an element of An-i and the homomorphism 

0n • An — >■ An-l 

is surjective and the kernel is The sequence 

... An^ An-i... -> A2^ Ai 

forms a "projective system" indexed by the integers > 1. Inverse 
limit of this inverse system is Zp = ^im(An, Refer [10] for details. 

Definition 1. The ring of p-adic integers, Zp, is the projective limit 
(inverse limit) of the system {An,4>n)- 

An element of Zp = lim(A„, 0„) is a sequence x = (..., x„, Xi) with 
Xn € An and (/>„(x„) = Xn-i if n > 2. Addition and multiplication in 
Zp are defined co-ordinate wise. In other words, Zp is a subring of 
the product nn>i ^n- If An is endowed with discrete topology and 
Y[n>i the product topology, the ring Zp inherits a topology which 
turns it into a compact space. 

Let p be a prime. For n G Z, let t'p(n) denote the exponent of 
highest power of p that divides n, and t'p(O) = oo. More formally, 
i/p(n) is the unique natural number such that n = p'^p^'^^u with 
p \ u. This definition can be extended to Q by letting t'p(f) = 
i'p{a) — i'p{h). The p-adic absolute value is defined by \x\p = p~^p^^'^ for 
any x G Q. Then dp{x, y) = \x—y\p defines a metric on Q. The metric 
space (Q, dp) is not complete, and its completion is the p-adic number 
field Qp. This absolute value is non- Archimedean as, in the place of 
triangle inequality, the stronger relation \x + y\p < Max{\x\p, \y\p}, 
also known as ultrametric inequality holds. The non-Archimedean 
property is equivalent to the assertion that, sup{|n| : n G Z} = 1 
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As a consequence (in stark contrast to the Euclidean norm), the 
p-adic norm does not permit accumulation of errors in the following 
sense: if each of k elements {xi,X2, ■■■,Xk} have p-adic norm atmost e, 
then \xi + X2 + ■■■ +Xk\p < e as well. This property justifies extensive 
use of modular arithmetic {p-adic estimation) in p-adic calculations. 

Definition 2. A p-adic integer is a formal series Y2i>o ^iP^ with 
integral coefficients ai satisfying < < p — 1. 

The subset of defined by Zp = {x G Qp : \x\p < 1} is called the 
integer ring and is an integral domain. The set of invertible elements 
in Zp, called its group of units is Z* = {a; e Zp : \x\p = 1}. The ring 
Zp contains a unique maximal ideal pZp = {x G Zp : |x|p < 1}. The 
quotient of Zp by this maximal ideal is a field, which can be identified 
with the finite field of p elements Fp = Z/pZ in the obvious way. Fp 
is called the residue class field of Qp. The field Qp is unorderable, 
in essence due to the modular arithmetic of Fp. It has characteristic 
since it contains Q as a subfield. Indeed, Q is a dense proper 
subset of Qp and embeds into Qp as the set of elements whose p- 
adic coefficients are eventually periodic. Topologically Qp is a Cantor 
set: totally disconnected but not discrete. In algebraic sense too p- 
adics are full of holes: There is no finite extension of Qp which is 
algebraically closed. Ostrowski's theorem states that any non-trivial 
absolute value on the rational numbers Q is equivalent to either the 
usual real absolute value or a p-adic absolute value. Refer [TT], [S] 
for details. 

3 p-adic Dynamics 

It is natural to have accumulated truncation errors or round-off er- 
rors in a dynamical system even with a small perturbation. These are 
unavoidable since even with simple repetitive operations, the number 
of digits of the result can increase so much that the result cannot be 
held fully in the registers available in the computer. Such errors accu- 
mulate one after another from iteration to iteration generating new 
errors. These difficulties motivated to look for an alternate number 
system which possesses the best features as well as the advantages of 
both the p-ary and residue number system. Such a number system is 
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the p-adic number system, discovered by Kurt Hensel in 1897 in the 
course of his work on finding new completions of the rational num- 
bers. Hensel's original description of the p-adic numbers involved an 
analogy between the ring of integers and the ring of polynomials 
over the complex numbers, the crux of which was the development 
of a representation of rational numbers analogous to that of Laurent 
expansions of rational functions namely, the p-adic expansion. This 
idea was motivated by the existence of real expansions of rational 
numbers with respect to a p-scale: 

k 

X = ^ «nP", a„ = 0, ...,p - 1. 

n=—oo 

Such manipulations with rational numbers and series generated the 
idea that there exists some algebraic structure similar to the system 
of real numbers M. Thus each has the structure of a number field. 
In fact, the fields of p-adic numbers, Q^, were the first examples of 
infinite fields that differs from Q, M, C and corresponding fields of 
rational functions. The following definition is useful in identifying a 
p-adic integer. 

Definition 3. Let p be some prime number. A sequence of integers 
satisfying 

Xn = a;„_i (mod p") 

for all n > 1, is called a p-adic integer. Two sequences and 
and {x'n} determine the same p-adic integer if and only if Xn = x'^ 
(mod p^) for all n > 

This definition can be easily identified with the definition [2l 

The connection between congruences and equations is based on the 
simple remark that, if the equation 

F{xi,X2, ...,Xn) = 

where F is a polynomial with integral coefficients, has a solution in 
integers, then the congruence 

F{xi,X2, ■■■,Xn) = (mod m) 
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is solvable for any value of the modulus m. On the other hand the 
situation is more complicated for congruences. For any modulus m > 
1, there are polynomial congruences having no solutions. For eg., the 
congruence — x + 1 =0 (mod m) has no solution if p is any prime 
factor of m by Fermat's theorem. Where as a congruence can have 
more solutions than its degree, for eg., — 7a; + 2 = (mod 10) 
has four solutions x = 3,4,8,9. But if the modulus is a prime, a 
congruence cannot have more solutions than its degree. 
The following results are basic, see [D]. 

Theorem 3. // the degree n of f{x) = (mod p) is greater than 
or equal to p, then either every integer is a solution of f{x) = 
(mod p) or there is a polynomial g{x) having integral coefficients, 
with leading coefficient 1, such that, g{x) = (mod p) is of degree 
less thanp and the solutions of g{x) = (mod p) are precisely those 
of f{x) =0 (modp). 

Theorem 4. The congruence f{x) = (mod p) of degree n has 
atmost n solutions. 

Corollary 1. // + + ... + feo = (mod p) has more 

than n solutions, then all the coefficients bj are divisible by p. 

Theorem 5. The congruence f{x) =0 (mod p) of degree n, with 
leading coefficient a„ = 1, has n solutions if and only if f{x) is a 
factor of x^ — X modulo p, that is, if and only if x^ — x = f{x)q{x) + 
ps{x), where q{x) and s{x) have integral coefficients, q{x) has degree 
p — n and leading coefficient 1, and where s{x) is a polynomial of 
degree less than n or s{x) is zero. 

The proofs of the above theorems are simple consequences of Fer- 
mat's little theorem and its application which can be found in many 
books on number theory. 

3.1 Hensel's Lifting Lemma 

As both Q and M are not algebraically closed, they do not always 
contain all roots of polynomials with integer coefficients. Though 
C is algebraically closed, as the degree of the polynomial increases, 
finding roots of the polynomial is computationally not feasible, as 
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it requires high level of precision for machine computation. To this 
end, we now turn our attention to solving polynomial congruences 
modulo prime powers. 

We note that for any polynomial f{x) e Z[x] and any integer r, there 
is a polynomial Qrix) G Z[x] with /(x + r) = f{r) + xf'{r) + x'^gr{x). 
This can be seen either through the Taylor expansion for f[x + r) or 
through the binomial theorem in the form 



The above standard results mentioned in the form of Theorems can 
be used to find the solutions to f{x) = (mod p). The question is 
how we might be able to lift a solution to one modulo for various 
exponents k. With regard to this, Hensel's Lemma is a powerful tool 
which relates the roots of a given polynomial to its solution modulo 
a prime. The lemma and its proof both rely on iterative procedures 
that return an agreeable solution if supplied with a well-behaved 
seed. 

Definition 4. If f{a) = (mod p), then the root a is called nonsingular 
if f'{a) ^ (mod p); otherwise it is singular. 

Two versions of Hensel's Lemma are stated below. 

Theorem 6. Hensel's Lemma over the ring of integers. 

Suppose that f{x) is a polynomial with integral coefficients. If f{a) = 
(mod pP) and f'{a) ^ (mod p), then there is a unique t (mod p) 
such that 

f{a + tp>) = (modp?+i). 

Theorem 7. Hensel's Lemma over the ring of p-adic integers. 
Let f G he monic. If oq E Tj is a simple root of f{x) = 

(mod p), then 3y E Zp such that y = (mod p) and f{y) = 0. 

Proof. Suppose that 3a„ such that /(a„) = (modp"). We must 
show that be lifted uniquely to a„+i (mod such that 

a„+i = a„ (mod p") and /(a„+i) = (mod p""*"^), then y is the limit 
of this sequence of (mod p'^) solutions. 




(1) 
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Since / is a polynomial we can write it in the form f{x) = Yli CiX^- 
Also consider tp" + possible lift of ■ Then 

/(a„ + tp") = ^Q(tp" + a„)^ (2) 

i 

= f{an)+p''tf\an) (modj9"+i). (3) 

The equivalence above is a result of Taylor series expansion. Now, 
solve for t in 

P"t/K) + /K)=0 (modp"+i). 

Thus 

t/'(an) = -(^) (modp). 

Since /(a„) = (mod p"), (since a„ = Qq (mod p)) and f\an) ^ 
(mod p) (simple root), then t has a unique solution (mod p). Thus 
a„+i = a„ + tp'^ is a unique lift of a„ (mod p). Thus we have con- 
structed an infinite sequence of such that /(aj) = (mod p*), such 
that, /'(flj) ^ (mod p*) and Oj+i = (modp). This sequence is 
Cauchy, and therefore converges to a unique limit y G Zp. 

Since Z is dense in Zp, proof of Theorem [6] follows directly from 
Theorem [71 

It can be verified that, proof of Hensel's Lemma is entirely analogous 
to Newton's method for locating the root of a differentiable function. 
Let us recall Newton's method from calculus as a method of finding 
roots to a polynomial by choosing a seed and then making better and 
better approximations based on the polynomial's derivative at that 
point. In the case of Newton's method, the condition on the seed is 
that the derivative at that point be non-zero, otherwise it supplies no 
useful information for improving at each iteration. Hensel's Lemma 
is similar, it takes a polynomial with coefficients in Zp and instead of 
requiring a guess at a possible root, it requires a p-adic integer that 
is a root mod p, i.e. some a such that the polynomial f{x) evaluated 
at a is 

/(a) = (mod pZp) 

This method will then return roots mod p,p'^,p^, ... until the desired 
root of the equation is found. 
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3.2 Backward Iterations, Inverse Limits and p-adic 
Approximations 

Now we find the backward iteration of any given polynomial f{x) of 
arbitrary degree say n, at any point, say Xn of the forward iterating 
orbit, i.e., we would like to solve 



for some A; G N. For this, we first solve the congruence f{x) = Xn 
(modp), such that the nonzero coefficients of f{x) are relatively 
prime to p and x„ is chosen as above, so that the corresponding 
backward iterating orbit of Xfi 8bTQ found modulo p^, j — )■ oo. 
If the degree n of f{x) is greater than p, then by Theorem |3l f{x) is 
divided by (x^ — x) mod p and solutions of the resulting polynomial 
g{x) are the same as those of f{x) = x„ (mod p). Since the modulus 
is prime, the congruence cannot have more solutions than its degree. 
Let ai,a2,...,ai be the roots obtained, by solving the congruence 
f[x) = Xn (mod p) where (/ < deg{g{x) by Theorem H]). Now each 
tti is lifted modulo p^ , j = 2,3,... whenever /'(oj) ^ (modp), 
{i = 1,2, ..,/.)i.e., whenever the a[s are nonsingular. We work up to 
a fixed precision say j = k. Once the roots are lifted modulo p^, Xn is 
replaced by the lifted root and the congruence equation (jl]) is soved 
now with respect to the new root. The process of replacing the old 
root by the new lifted root is repeated for some finite number of steps. 
Thus with respect to each nonsingular a^, we obtain the correspond- 
ing backward sequence, generated from the single seed Thus if 
there are n nonsingular roots, after the mth step of replacing the old 
root by the new lifted root, there are atmost backward iterating 
points generated from the single seed Xn- Hence there exists a tree 
like structure, the roots may be called as leaves and the branches are 
formed at each new lifted root. A simple code for this program can 
be written on python, which computes the sequences effectively. 
Sequences thus obtained belong to the inverse limit space by defini- 
tion. It can be verified that the sequence space formed by the above 
backward iterations is totally disconnected and discrete. Also, if the 
original spaces are discrete, then the inverse limit space ^m{A, /} is 
totally disconnected. This is one way of realizing the p-adic numbers 
and the Cantor set. 




(4) 
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To study the long time behaviour of a dynamical system it is nec- 
essary to introduce a suitable metric. A natural choice would be 
the one given in section [H i.e., the distance between two sequences 
s = (so, si, ...) and t = (to, ...) is given by d{s, t) = X;,=o 
Thus by the introduction of a metric, sequence space of backward 
iteration becomes a compact metric space. The sequence space thus 
obtained can be identified with the ring of p-adic integers in view of 
the following theorems: 

Theorem 8. Any two totally disconnected perfect compact metric 
spaces are homeomorphic. 

Theorem 9. Let M be a compact , totally disconnected metric space. 
Then M is homeomorphic to the inverse limit space of an inverse 
limit sequence of finite discrete spaces JBi- 

4 Conclusion 

The problem of backward iteration is studied by solving congruences, 
which in turn are solved by Hensel's lifting. Sequences generated by 
such solutions form, naturally, elements of an Inverse Limit Space. 
These spaces have also been characterised. 
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